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Running coupling constant and propagators in SU{2) Landau gauge* 

Jacques C. R. Bloch^, Attilio Cucchieri^ Kurt Langfcld'' and Tereza Mendes''^ 
'^Institut fiir Theoretische Physik, Universitat Tubingen, D- 72076 Tiibingen, Germany 
''IFSG Sao Paulo University, CP. 369 CEP 13560-970, Sao Garlos (SP), Brazil 

We present a numerical study of the running coupling constant and of the gluon and ghost propagators in 
minimal Landau gauge. Simulations are done in pure SU{2) lattice gauge theory for several values of /? and 
lattice sizes. We use two different lattice setups. 



1. INTRODUCTION 



We consider, on the lattice, a running coupling 
constant g^{p) defined by ||,|^ 



Hp) 



9l 



(1) 



where D{p) and G{p) are, respectively, the gluon 
and ghost propagators evaluated in Landau 
gauge. Clearly g'^{p) is a gauge-dependent quan- 
tity; however, notice that g^ip) is renormaliza- 
tion-group invariant in Landau gauge since, in 
this case, ZgZ^ Z3 = Zi = 1. This running cou- 
pling strength enters the quark Dyson-Schwinger 
equation directly and can be interpreted as an 
effective interaction strength between quarks |^ . 

Studies of the coupled set of Dyson-Schwinger 
equations for the gluon and ghost propagators 
have shown that: (i) the gluon propagator be- 
haves as D{p) ^ p-2+4K jj-^ -(.jjg infrared limit 
[and thus D{0) = if k > 0.5], (ii) the ghost 
propagator behaves as G{p) ~ |,-2-2k small 
momenta and (iii) the running coupling strength 
ttsip) = g^(p)/47r defined in eq. ([l]) has a finite 
value ac at zero momentum (infrared fixed point). 
Using different approximations, in order to solve 
the Dyson-Schwinger equations, the following val- 
ues have been obtained: k « 0.92 and ac ~ 9.5 
K w 0.77 and « 11.5 @, k w 0.60 and 
ac ~ 8.9 /Nc 0. [Here, the first two results re- 
fer to SU{3).] We stress that the large value 
for ac obtained in [Qj^ is related to the angu- 
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lar approximation used in the integration kernels. 
Let us notice that, using stochastic quantization 
||^, Zwanziger also obtained that the transverse 
gluon propagator in the infrared limit behaves as 
Dip) p-2+4«: ^^^Yl K w 0.52. 

From the lattice point of view we know that lat- 
tice gauge-fixed Landau configurations belong to 
the region delimited by the first Gribov hori- 
zon, and that Q is not free of Gribov copies. One 
can also prove [^ that the restriction of the path 
integral to the region implies a suppression of 
the (unrenormalized) transverse gluon propaga- 
tor D(jj) in the infrared limit. At the same time, 
the Euclidean probability gets concentrated near 
the Gribov horizon and this implies enhancement 
of G{p) at small momenta m. 



2. RESULTS 

Simulations have been done in Sao Carlos for 
/? = 2.2, 2.3,..., 2.8 and V = 14*, 20^ 26^ 
and in Tubingen for f3 — 2.1, 2.15, 2.5 and 
y = 12^ X 24, 16^ X 32. The simulations carried 
out in Tiibingen are based on a direct evaluation 
of the form factors F{p) = D{p)p^ and G{p)p^ 
appearing in eq. (0). Also, for the evaluation of 
F{p), the gluon field has been defined in terms 
of the adjoint links [|| instead of the usual link 
variables. The gluon field obtained in this way is 
invariant under non-trivial Z2 transformations. 

Gribov-copy effects for the two propagators, if 
present, are smaller than the numerical accuracy 
9|. Preliminary results have been presented in 

io|. 
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Figure 1. Fit for the running coupling using eq. 
(|) with Co = 1.4(2), aa ^ 5.5(3), S = 1.77(9), 
A = 0.83(4) and A set to 2.2. 
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Figure 2. Fit for the ghost and gluon propagator 
form factors using eqs. (||) and (||) respectively, 
with ci = 0.98(4), C2 = -0.59(6), A = 0.98(2), 
B — 1.124(9) and a{p) as obtained from the fit 
reported in Fig. Ill 



In order to compare lattice data obtained for 
the two propagators at different (3 values we 
used a standard scaling analysis ||ll| based on 
maximum overlap without considering any phe- 
nomenological fit functions. (Details will be pre- 
sented in [|l2j .) Also, for the data produced in Sao 
Carlos, we have discarded data points at small 
momenta that are affected by finite-size effects. 
(These finite-size effects are less pronounced when 
one evaluates the form factor directly.) 

We have considered two different sets of fitting 
functions, namely 



a{p) 



1 



co + t' 



D{p)p^ = A- 



coflQ -I- a2{t + A) 
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13/22 



Cl 



G{p)p' 
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Cl 
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(2) 
(3) 



t 



,'/'\p) (4) 



where t = p^/A^ and 02 (p) is the 2-loop running 
coupling constant IQ, and 



a{p) ^ Cp*/ [(/ + m) s{a)] 
D[p) = Ap' / [{p^ + m) s-<" [an)] 
G{p) = Bj {p^s^^{aG)\ 



(5) 
(6) 
(7) 



where s(a) = (ll/247r2) log [1 + (pVA^)"], 7^ = 
13/22 and 70 = 9/44. Note that, in the first case. 



the fitting functions correspond to k — 0.5, while 
in the second case one has kq = aclG and kd = 
1 — aDlDl'i'- Also, both sets of fitting functions 
satisfy the leading ultraviolet behavior of the two 
propagators. 

Results of the fits are reported^ in Figs. ||-||. 
From our data there is evidence for the suppres- 
sion of the transverse gluon propagator D{p) in 
the infrared limit and for the enhancement of the 
ghost propagator G(j)) in the same limit. Also, 
the running coupling strength as(p) defined in 
eq. (H) probably has a finite value at zero mo- 
mentum. However, in order to probe the infrared 
region and give a final value for k and olc one 
needs to simulate at larger lattice volumes. 
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Figure 4. Fit for the gluon propagator using eq. 
d) with A = 1.02(9), ao = 1.9(3) and m = 
0.8(3); this gives kd ^ I - aDlo/'^ ^ 0.44(9). 
Here A has been set to 1.322 (see Fig. ^). If 7d 
is also a fitting parameter we get 7_d — 0.579(7) 
to be compared with 13/22 « 0.591. 
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Figure 5. Fit for the running coupling a{p) us- 
ing eq. (|) with C = 0.072(8), a ^ 1.9(3), 
A = 1.31(1) and m = 1.0(6). 



